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Abstract 

A massless up quark is an intriguing possible solution to the strong 
CP problem. We discuss how lattice computations can be used in 
conjunction with chiral perturbation theory to address the consistency 
of m u = with the observed hadron spectrum and interactions. It 
is not necessary to simulate very light quarks — three flavor partially 
quenched computations with comparable sea and valence quark masses 
on the order of the strange quark mass could suffice. 
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1 Introduction 



QCD allows violation of the symmetry CP through the parameter 9 = 
Qqcd — arg det.M, where M. is the quark mass matrix. All experimental 
evidence implies that 6 < 10~ 9 , and thus CP is nearly preserved by the 
strong interactions. The absence of a symmetry to ensure 9 = is one of the 
most perplexing features of the Standard Model. Three different solutions 
are commonly considered for this "strong CP puzzle" : dynamical relaxation 
of 9 by means of an axion field [|IJ; spontaneous CP violation at high energy 
with a Nelson-Barr mechanism 0, |3] to ensure reality of the quark mass de- 
terminant; and a vanishing up quark mass m u = 0. In the particularly simple 
case of vanishing m u the determinant of the quark mass matrix is zero, and 
9 is no longer a physical parameter. Models in which the up-type quark 
mass matrix has rank two at short distance as an accidental consequence of 
symmetry are easily constructed @, H H 0, |[- 

The viability of a massless up quark can in principle be determined by 
comparing the predictions of QCD with the observed spectrum and interac- 
tions of hadrons. However, effective hadronic theories whose sole input from 
QCD is the approximate SU(3) x SU(3) chiral flavor symmetry cannot settle 
this question. Under SU(3)l x SU(3)r, the quark mass matrix M. trans- 
forms as the (3, 3) representation, and therefore so does the matrix 1 1 M. \ . 
Thus, symmetry considerations alone can never rule out the possibility that 
explicit symmetry breaking occurs in hadronic physics through combinations 
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of the form M. e ff = M.+M.^ 1 |7W | /A, with A a scale determined by strong 
interactions, an ambiguity pointed out by Kaplan and Manohar ||. In the 
case m u = 0, this gives an effective quark mass matrix 



M 



eff 



I m* d m*/A \ 
m d 

V m 8 ) 



(1) 



which allows the quantity m d m*/A to fulfill the role conventionally played by 
m u . Note, however, that in this case det M. e ff is rea l> an d there is no strong 
CP problem. The nonlinear term in M. e ff could arise from instantons, for 
example, where the d and s quark zero-mode propagators are connected by 



m d and m s insertions respectively P, 10, 11]. Appropriate values for the A in 



this scenario would correctly fit all hadron data. The conventional extraction 



of quark mass ratios from chiral perturbation theory would nevertheless 
incorrectly yield a nonzero value for m u /m d . 

To resolve the ambiguity chiral symmetry may be supplemented with ad- 
ditional assumptions. In Ref . |H| , Leutwyler gives a list of three plausible 



assumptions, each of which independently rules out m u = 0. The assump- 
tions are: SU(3) symmetry is always approximately valid, and all physical 
quantities can be reliably expanded in powers of m s ; dispersion relations are 
saturated by the lowest lying states; and the large- A^ c explanation for Zweig's 
rule (suppression of virtual quark loops by 1/AQ is valid. 

We consider the above arguments against m u = to be quite reasonable. 
Each of these assumptions can be experimentally tested in various ways, and 
none has yet been disproven. Yet it is possible that an effective up quark 



mass, although sub-leading in 1/N C , is numerically large enough to account 
for the hadron spectrum without implying a general breakdown of chiral 
symmetry. Due to the importance of the strong CP problem, we consider it 
essential to pin down the value of m u without assumptions. 

Lattice simulation of QCD can, in principle, decide whether or not a 
massive up quark is required, particularly in light of recent advances in re- 
alizing chiral symmetry on the lattice [|14| , [15], [16|. With the aid of chiral 



perturbation theory, it is not necessary to simulate QCD with a massless or 
extremely light quark. We begin by reviewing chiral perturbation theory in 
the continuum, and then we consider in turn full, partially quenched, and 
quenched QCD on the lattice. We find that lattice simulations in either full 
QCD or partially quenched QCD with all quark masses comparable to the 
strange quark mass could largely settle the issue. 

2 m u in the continuum 

The light quark masses (m u , m^, m s ) are small compared with the character- 
istic mass scale of the strong interactions, and QCD possesses an approximate 
SU(3) CED SU(3) chiral symmetry broken spontaneously to the vector SU(3) 
subgroup, resulting in a light pseudo-Goldstone boson octet with predictable 
low energy interactions [0, [17], [L8|, [L9|, [2C], [2l| . For sufficiently light quarks, 
the pseudo-Goldstone masses squared are nearly linear in the quark masses, 
and the quark mass ratios may be extracted from the pseudoscalar spectrum 
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with values m u /m d = 0.56, m s /m d = 20.1 |T7]5 However, the strange quark 
mass is large enough that the lowest order predictions from chiral symmetry 
receive significant corrections, of size m s /A x , where A x is the chiral symme- 
try breaking scale of order a GeV. The simplest way of extracting predictions 
from chiral symmetry to any given order in m s and in pion momenta is to 
use a phenomenological chiral Lagrangian Using the parametrization 



of Gasser and Leutwyler (JT2J, [L^, |2C| ^TJ the Lagrangian relevant for the 
extraction of quark masses to second order in m s is 



C 



c 2 + u + ■ ■ ■ 

£ Tr (d^u) + £ Tr (^U + X U ] ) 
. . . + U Tr (dprfdPTj) Tr {u ] x + X ] u) 
+L 5 Tr IdM^d^U (U ] x + X ] u) 



Tr (rt x - X ] U 



+L,[Ti(u^x + x ] U)] 2 + L, 
+L 8 Tr (x ] U X ] U + xU ] xU^) + . . . + L 12 Tr X ^ 



(2) 
(3) 



(4) 



where 



M 

U 

X 



diag(m U} m d ,m s 
exp 



i2T a ir a 



f 



= 2MB 



(5) 
(6) 
(7) 



T a are SU(3) generators, n a are pseudoscalar fields, and B,L iy f are phe- 
nomenological parameters characterizing QCD dynamics. (It is not possible 
1 Similar earlier quark mass estimates were given in Ref. [E2|. 



to separately determine M. and B from experiment.) Most of these parame- 
ters are renormalization scale dependent — we follow the usual practice in the 
chiral perturbation theory literature of quoting all parameters at the p mass. 
They are then easily run to a different scale by using the renormalization 
group. For convenience electromagnetism has been left out of this effective 
theory although electromagnetic effects contribute to the pseudoscalar me- 
son masses. We therefore define "QCD" masses which have electromagnetic 
effects subtracted — to leading order these are just the physical meson masses 
with the exception of the electrically charged mesons |23| 



Ml ±QCD « il4 2 0phvs (8) 



L ir±QCD ~ JW 7r°phys 

r2 r ^ J T\/r2 i\/r2 , tij-2 



M 2 K±QCD « M 2 K±phys - il^ ±phys + M^ 0phys . (9) 

These lowest order formula for the QCD masses can be improved beyond 
leading order [^4), |25], |28], |29|. All meson masses used in subsequent 



formulae are these QCD masses. From this effective theory, it is possible to 
determine a constraint on quark mass ratios, up to corrections of order m 2 s . 



20, 21, 30 



m * ~ ™° =02^2 = M% {Ml -Ml) 

mf - ~ Q ~ Qd ~ Ml (M| - M| ± ) ~ (242) ' (10) 

where 

rh= \ (m u + m d ) . (11) 

If the values of the Us are known, the quark mass ratios may be found: 

m s 2Ml , , 

A 1 (12) 



m M2(l + A 



M) 
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where 



m„ 



m 



Ml M 2 K - (1 + A M )M 2 

g 2 (i + A M ) 2 



(13) 



A M = -fr + n n + —{Ml - M 2 )(2L 8 - L 5 



(14) 



and 



M 2 



32tt 2 / 2 



log 



'Ml 



(15) 



To decide whether m u can be zero, the value of 2L 8 — L 5 is required. L 5 can be 
determined from the ratio fx/fn, with a value 10 3 L 5 w 1.4 ±0.5. The linear 
combination L5 — 12Lj — 6L§ can be extracted from the pseudoscalar masses, 
giving 10 3 (2L 7 + L 8 ) = ±0.2. (The quoted errors represent the estimated 
theoretical uncertainty due to higher order corrections in m s .) In principle 
(but not in practice) L 4 and L 6 — L7 can be determined from meson-meson 
interactions. However data cannot be used even in principle to completely 
specify the L's. The Gasser-Leutwyler Lagrangian is invariant under the 
replacement |21] 



Lq — 6 Lj — > L7 — 5 L§ — > L§ + 26 



(16) 
(17) 



with 



6 = /3/ 2 /32 . 



(18) 



This invariance is a necessary consequence of the Kaplan-Manohar ambiguity, 
and leads to a corresponding ambiguity in the L's. 
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A combination of phenomenological and large N c constraints can be used 
to determine all the L's with an estimated 30% precision [IT], The 
results agree (by design) with the picture that all terms in £ 4 are consistent 
with the values obtained by integrating out light resonances, and are incon- 
sistent with m u = . In particular, a conventional estimate of this kind 
gives H, H, H, |6| 10 3 L 8 « 0.9 ± 0.3, 10 3 L 7 « (-0.4 ± 0.2), while the 



hypothesis m u = would require 10 3 L 8 = -0.4 ±0.3, and 10 3 L 7 = 0.2 ±0.2. 
Thus any method which can directly determine the L's, even with errors as 
large as 100%, can distinguish between these two possibilities. 

3 m u on the lattice 

Simulating QCD on a lattice is, in principle, the most reliable way of de- 
termining quantities which are not predicted from symmetry alone. For in- 
stance, although it is only possible to constrain ratios of light quark masses 
phenomenologically, recent lattice computations of the hadron spectrum have 



suggested an absolute range for the strange quark mass |37|, [38|, f4(| |41] 



42[ [43| , |44] , |45| , f46| and the parameter B. It is difficult to simulate quarks 
with realistic masses: light quarks require large lattices to avoid finite vol- 
ume artifacts. However with the aid of chiral symmetry lattice computations 
done at moderate quark masses can be extrapolated to lighter quark masses. 
Unfortunately this extrapolation is complicated and plagued with unphysi- 
cal artifacts in the usual quenched approximation, where quark dynamics are 
not included. 
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3.1 Unquenched lattice QCD 

With computing power sufficient to simulate unquenched QCD with quark 
masses light enough to apply chiral symmetry, there are several possible 
methods^ to extract the chiral Lagrangian coefficients and determine m u . 
One way is to use the chiral symmetry prediction for the pion mass squared 



Ml 



2mB ( 1 + - + 2mK 3 + -± ^ S -K 4 



(19) 



where 



^3 = 

K 4 = 



8B 

1 
485 



2 (2^8 - L 5 ) 



/ 2 



(2L 6 - L 4 ) 



(20) 
(21) 



By measuring the pion mass as a function of rh and m s the combina- 
tion 2Lg — L5, needed in determining the quark mass ratios, as well as the 
combination 2L$ — L4, which provides an interesting test of conventional 
assumptions and the large N c expansion, may be extracted. Independently 
varying rh and m s requires a lot of different simulations, however. It may be 
simpler to work with equal quark masses m u = rrid = m s = rh = rh and vary 
rh. A fit to the quadratic dependence of the pion mass squared then yields 
the linear combination 2L$ — L5 + 6L$ — 3L4 from 



d 2 Ml 
dm? 



127T 2 / 2 



3 + 2 log ( + 768tt 2 (2L 8 - L 5 + QL 6 - 3L 4 ) 



V 



.(22) 



'For earlier work on lattice computations of chiral coefficients see rcfs. HI 
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With equal quark masses, the combination 2L 6 — L 4 may be separately ex- 
tracted by a measurement of the matrix element of ss in the pion, using 



.... mB 2 



36tt 2 / 2 



l + logf^Ull527r 2 (L 4 -2L 6 ) 



(23) 



V J 

Lq may be independently extracted from the vacuum expectation value (/ ss J dd). 
and Li may also be measured from the dependence of /„- on fh. Thus in prin- 
ciple by simulating QCD on the lattice with several different quark masses it 
is possible to verify the conventional estimates of those L's for which there 
is no direct experimental data. Even a measurement of 2L 8 — L 5 with 100% 
errors provides an interesting test of the large N c expansion and, depending 
on the result, could rule out the possibility that m u = 0. Note that mea- 
surement of Li and Lq would also be quite interesting theoretically, although 
not directly needed to extract quark mass ratios. Available data provides no 
constraints on L^q. 

Ruling out m u = along the above lines may not be easy. Here we define 
the more restrictive "effective up mass hypothesis," which may be somewhat 
simpler to test than whether m u = 0, since this hypothesis makes a prediction 
for L 6 . The effective up mass hypothesis is motivated by the agreement 



between different experimental determinations of light quark mass ratios (T2|, 
|3C| p9| , |50| , |51~1 . The hypothesis is that m u = 0, but that the conventional low 
energy theory works accurately with the replacement m u —>■ m* d m*/h. In the 
chiral lagrangian, such an effective up mass is a nonstandard contribution to 
the coefficients Lq,Lj, and Lg, in the combination ALq=ALj = —2AL$ = 
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0.7 x 10~ 3 . With this hypothesis, all the L's are determined, and L 4 agrees 
with the conventional estimate 10 3 L4 = —0.3 ± 0.5, while 10 3 Lg is 0.5 ± 0.3, 
in contrast to the conventional estimate 10 3 L 6 = —0.2 ± 0.3. 

3.2 Partially quenched lattice QCD 

QCD simulations in the "partially quenched" approximation, which includes 
the effects of N flavors of dynamical quarks with mass m sea different from 
the valence quark mass, may also be of interest. Such simulations could 
be decisive for the determination of the Gasser-Leutwyler L coefficients 3 . 
For such an analysis to be reliable it is necessary that the sea quark mass 
is sufficiently small for the dominant artifacts of partial quenching to be 
computable using partially quenched chiral perturbation theory f5"4| . This 
requires that the mass M$s of a pion made of sea quarks be light compared 
with the scale A x . It is also desirable to take the valence quark mass to be 
comparable to the sea quark mass, to reduce quenching artifacts from the 
non-decoupling of the rj' |55|]0- In the limit m sea ~ m « A x , Sharpe has 
calculated the following dependence of the pion masses on the valence and 



sea quark masses |54 



Ml ± = 2mB{\ + — — [(2m - m sea ) log {^^-) + m - m sea 

L 1\87T Z J Z 1 jl l 

-(2L 8 - L 5 ) + t^(2L 6 - L 4 )} (24) 



,2mB s 



3 After the completion of this work we were informed of the work of Sharpe and 
Shoresh |32| who reach similar conclusions. 

4 Note added in revision: Recent work by Sharpe and Shoresh [^3| has shown that the 
artifacts from the rj' are under theoretical control even when the valence quark mass is 
much lighter than the sea quark mass, provided both masses are sufficiently small. 
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M SV = ( m sea + m)B 



(2mB\ 8(m sea + m)B MT 

1 + AWT log ) + p (2Ls " Ls) 



+ ^^(2L 6 -L 4 )1 (25) 



M 2 SS = 2m se& B 



m SP ,B , (2m RPa B s 



1 + »ro 9X9 lQ g 



N8ir 2 f 2 °\ /j, 2 



" (2L 8 -L 5 ) + ^£^(2L 6 -L 4 ) . (26) 



Here all valence quarks have mass fh, 7r ± is a pion made of different valence 
quarks, M$s is the mass of a pion made of sea quarks and M$v is the mass 
of a pion made of one sea and one valence quark. For N = 3 the parameters 
B and Li are the same as those in the QCD chiral lagrangian [pi, 56]. Thus 
for N = 3, the desired combination 2L§ — may be extracted by fitting the 
pion masses as a function of fh with m sea held fixed. Note that an iV = 2 
simulation, while interesting, is not sufficient to determine the L coefficients, 
as these may have significant dependence on the number of flavors. For 
instance the contributions from gauge field configurations with fermion zero 
modes, such as instantons, should be quite sensitive to the number of sea 
flavors. 



Golterman and Leung |55| have extended the partially quenched chiral 
perturbation theory calculations to the case where the ?/ is light compared to 
the scale A x , as expected in the large N c limit. In this limit, unless the valence 
and sea quark pion masses are comparable and both much lighter than the 
7]' , the pion masses depend on two new parameters associated with the rf 
mass and decay constant. Even with a light rf it is theoretically possible, 
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with enough different measurements, to extract 2L 8 — L 5 from lattice data[]. 
3.3 Quenched lattice QCD 

The quenched approximation is not a systematic approximation to QCD. 
Nevertheless it is generally used to facilitate lattice computations with cur- 
rently available computing power, and quenched QCD does have a spectrum 
similar to the real thing ^jTl [58]. This suggests that the dominant effects of 
quark loops can be compensated for in the quenched approximation by ad- 
justing the QCD parameters (quark masses and the scale Aqcd)- However 
a reliable extraction of the true value of the chiral Lagrangian coefficients 
and of m u from the quenched spectrum is problematic; the possibility that 
the quenched approximation with a nonzero up quark mass mimics the true 
QCD spectrum with m u = cannot be ruled out. For instance the quenched 
approximation is missing the down and strange quark loop effects which, in 
conjunction with instantons, might mimic an effective up quark mass pro- 
portional to m* d m*. 

A possible way to explore the quark loop contribution to the effective up 
quark mass in the quenched approximation is to explicitly include sources 
for the sea quarks. For instance, one could measure a three-point function 

(tt| Jss\n) . (27) 

In particular, the instanton effects which might give an effective up quark 

5 Note added in revision: see recent work of Sharpe and Shoresh |5^] for how to deal 
with the r( artifacts. 
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mass do contribute to this matrix element in the quenched approximation. 
In full QCD this matrix element is equivalent to dM%/dm s , which would 
vanish ignoring quark loops; however the equivalence does not hold in the 
quenched approximation. Logarithmically enhanced quenched artifacts from 
an rj' loop0 give a contribution to this matrix element which are suppressed 
by 1/iV? but which introduce a new parameter which cannot be computed 



using quenched chiral perturbation theory |59|, pOfl . To avoid this artifact it 
may be better to measure 

(tt| J {2ss - dd-uu)\n) . (28) 

4 Summary 

Chiral perturbation theory makes possible the use of lattice simulations of 

full and partially quenched QCD with moderate quark masses to learn about 

the properties of QCD with light quarks. In this paper we showed how to 

use full or partially quenched simulations with equal, moderately sized quark 

masses to extract the second order coefficients in the pion chiral Lagrangian. 

Such calculations are of interest to check the predictions of large N c QCD, 

to verify chiral perturbation theory, and to test the hypothesis of resonance 

saturation of dispersion relations. Such computations, even with large errors, 

can provide a method for settling the important issue of whether the m u = 

solution to the strong CP problem is consistent with the spectrum of light 

pseudoscalar mesons. 

6 We thank Steve Sharpe for explaining this to us. 
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